Electron Spectra for D and B Meson Semi-leptonic Decays at RHIC from PYTHIA 

with Modified Heavy Quark Fragmentation 
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We present a study of charmed hadron transverse momentum distribution from p+p collisions at 
yi]^ = 200 GeV using the PYTHIA Monte Carlo event generator. The PHENIX tuned PYTHIA 
parameter set based on electron measurements alone yields a D meson pr distribution softer than 
the preliminary STAR experimental data from d+Au collisions scaled by the number of binary 
collisions. In order to match the STAR pr spectrum, the higher order QCD effects in PYTHIA 
have to be taken into account, and a fragmentation function of charm quarks much harder than the 
default Peterson function is needed. Electrons from beauty quark semi-leptonic decays are found 
significantly below the contribution from charm quark decays for the pr region up to 8 GeV/c in our 
modified fragmentation scheme, while in the default fragmentation scheme the B decay electrons 
dominate at pr above 4 GeV/c. We propose an experimental measurement of electron and charged 
hadron correlation to differentiate these scenarios. 

PACS numbers; 24.10.Lx, 14.40.Lb, 13.20.Fc 



Heavy quarks are believed to be produced through 
initial parton-parton, mostly gluon-gluon, scatterings in 
nuclear collisions at the Relativistic Heavy Ion Collider 
(RHIC) energies. Theoretical calculations of heavy quark 
production within the perturbative Quantum Chromo- 
Dynamics (pQCD) framework is considered more reliable 
because the heavy quark mass sets a natural scale for 
the pQCD. Charm quark production is sensitive to the 
incoming parton flux for the initial conditions of nuclear 
collisions The transport dynamics of the heavy 

quarks in nuclear medium such as flow |3| and energy 
loss 00 can probe QCD properties of the dense matter 
created in nucleus-nucleus collisions. Therefore, heavy 
quark measurements provide unique insights into QCD 
properties of the new state of matter produced in nucleus- 
nucleus collisions. 

The heavy quark production in p-l-p or p-l-A collisions 
provides a reference for heavy meson formation and for 
nuclear modification factors of heavy quarks in nuclear 
medium. The charm fragmentation function has been 
measured in e^e^ and jp collisions and was fit by the 
Peterson function 6] (defined later in the paper) with a 
parameter e « 0.05 (default in PYTHIA 0) M- How- 
ever, in charm hadroproduction, it was observed that the 
c-quark px distributions of next-to-leading-order (NLO) 
pQCD calculations agree well with the measured open 
charm pT spectrum indicating that a much harder 
fragmentation function peaked at z « 1 is needed in 
charm hadroproduction. A more detailed discussion of 
this observation can be found in Coalesence 0| or 
recombination 0, models have also been proposed 
for charmed meson formation by combining a charm 
quark with a light up or down quark, presumably with 
soft pt [l^ . Thus the charmed hadron px would coin- 
cide with the bare charm quark px distribution in this 
hadronization scheme. These various hadron formation 
schemes can lead to significantly different interpretation 



of electron px spectra from experimental measurements. 

Recently the PHENIX collaboration reported a mea- 
surement of non-photonic electrons from heavy quark 
semi-leptonic decays from pp and A A collisions ^3 ■ 
The measured electron yields appear to follow approx- 
imately a binary scaling. Such a scaling would imply 
that the charm quarks have a small energy loss in the 
dense medium created at RHIC. However, such an inter- 
pretation is subject to the ambiguity that the electron 
contribution from B quarks, which is predicted to have a 
much smaller energy loss than that of charm quarks |0| , 
is unknown experimentally in the px region above 3-4 
GeV/c. In order to evaluate heavy quark energy loss 
using non-photonic electron measurement, contributions 
from D and B meson decays have to be investigated. 

In this paper we evaluate the px distribution of 
D mesons from PYTHIA v6.22 and compare the 
PYTHIA results with the STAR preliminary measure- 
ments p^. The charm quark fragmentation function 
will be modified from the default Peterson function and 
the PYTHIA parameters are tuned in order to match 
the PYTHIA results on px distribution with the exper- 
imental data. The electrons from D meson decays are 
compared with STAR preliminary electron px distribu- 
tion We present a method to qualitatively estimate 
the significance of B meson decays to the electron spec- 
trum at high Px and to verify the scheme of modified 
heavy quark fragmentation function as an explanation 
for hard px distribution of D meson in pp collisions at 
RHIC. 

In our PYTHIA calculations, we started with a 
PHENIX tuned parameter set 0| , and changed the value 
of PARP(67) to 4, which enhances c-quark production 
probability through gluon splitting ^20.] . In addition, 
we also modify the parameter in the Peterson function 
so that PYTHIA calculation approximately follows the 
shape of the STAR measured open charm px spectrum. 
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FIG. 1: and c-quark spectra from PYTHIA calculations 
compared with STAR preliminary data from d+Au collisions 
scaled by Nbin = T.5. 



Fig. n shows our PYTHIA calculations together with 
the STAR open charm spectrum [18| . where the mea- 
sured open charm data points from d-|-Au collisions at 
^/SNN = 200 GeV have been divided by Nbin — 7.5 - 
the number of binary collisions. The PYTHIA spectra 
have been normalized to the measured dN/dy for I?" at 
mid-rapidity (0.0265 ± 0.0036(stai.)±0.0071(s?/st)). 

The charm quark px distribution labeled by stars in 
Fig. mis from PYTHIA calculation with PHENIX tuned 
parameters and PARP(67)=4. We will refer this set of 
parameters as high-py-tuned parameters in the rest of 
the paper. The PARP(67)=4 was introduced to take into 
account higher order effects in pQCD calculation. We 
found that the change mainly affects c-quark production 
at high pt- Due to the large uncertainties in scaling D* 
and D± to D° {D*/D° « = 0.4±0.09±0.13), the 

STAR measurement of D pr shape has large systematic 
uncertainties in the high pT region. Without fine tuning 
any other PYTHIA parameters, we found that the bare 
charm quark spectrum approximately match the STAR 
open charm px distribution. This also indicates that 
PARP(67)=4 somehow does generate a pT distrbution 
similar to the actual NLO pQCD calculation .21J . which 
coincides with the STAR pT spectrum as well [ig . It is 
not the purpose of this paper to show that the PYTHIA 
calculation with our parameter set and NLO pQCD cal- 
culation are equivalent. The PYTHIA calculation with 
the high-pT-tuned parameter set seems to reproduce the 
same observation seen when comparing NLO pQCD cal- 



FIG. 2: (a) Electron spectra from PYTHIA calculation with 
the high-pT-tuned parameters and modified Peterson func- 
tion for charm quarks and beauty quarks compared with 
background subtracted single electron spectrum measured by 
STAR EMC from d-f-Au collisions scaled by Nbi„ = 7.5. (b) 
Electron spectra from PYTHIA calculation with the high-pT- 
tuned parameters and default Peterson function for charm 
quarks and beauty quarks. 

culation for c-quark to measured open charm data [lOj . 
We will further modify the fragmentation function in PY- 
HTIA and investigate the consequences on single electron 
measurements . 

The solid squares in Fig. ^ depict pT distribution 
from PYTHIA calculation with the high-px-tuned pa- 
rameters and the default Peterson fragmentation func- 
tion in PYTHIA. The default Peterson function refers to 
the value of the parameter e in Peterson function 

D{z=pDlPc)c^-r —. ^—rj- ^ (1) 

z|l — 1/z — e/yl — z)Y 

being 0.05 for charm quarks and 0.005 for beauty quarks. 
The default Peterson fragmentation for charm quarks is 
too soft to reproduce the measured open charm spec- 
trum. We modified the value of the parameter e to make 
the fragmentation function harder. The result is shown 
as open circles in Fig. ^ The Z?" spectrum roughly 
coincides with the measured data. Here the value of 
£ we used is 10~^ for both charm quarks and beauty 
quarks. In this case the fragmentation function is very 
like (5(1 — z). Hence, the charm quark hardly loses its 
momentum during fragmentation. The PYTHIA calcu- 
lation with high-pT-tuned parameter set and modified 
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FIG. 3: Sum px distributions for charged hadrons in the 
electron cone. Solid lines show electrons from B meson decay 
case. Dashed lines show electrons from D meson decay case. 
Left panel: Electrons selected at pr > 2.0 GeV/c. Middle 
panel: Electrons selected at pr > 3.5 GeV/c. Right panel: 
Electrons selected at pr > 5.0 GeV/c. 



Peterson fragmentation function (e — 10^^) can reason- 
ably depict the measured open charm data. A harder 
fragmentation function is needed for the hadronization 
of charm quarks if the pQCD calculation is to describe 
the STAR preliminary data. 

The STAR independent measurements of the recon- 
structed D'^ and single electrons from heavy quark 
semileptonic decays measured with TOF and TPC are 
consistent The electron measurement there only 

covers up to pt < 4 GeV/c and has no sensitivity to 
the B meson contribution. We also checked the consis- 
tency between open charm data and non-photonic sin- 
gle electron data in d-fAu collisions measured by STAR 
EMC [13 within our PYTHIA calculation. The results 
are shown in Fig. |21 The PYTHIA spectra of electrons 
from charm meson decays are scaled by the same fac- 
tor used to scale the PYTHIA D'^ spectra to the mea- 
sured dN I dy{D^) at mid-rapidity. The electron spectra 
from beauty meson decay are normalized by the ratio of 
^bbl^cc based on the NLO pQCD calculation The 
band corresponds to the theoretical uncertaint y o f this 
ratio, which is in the region of 0.45% - 0.60% y^. We 
used the value at the center of this range to calculate the 
sum of the electrons from both charm and beauty me- 
son decays. The EMC electron data points are divided 
by Nbin- The electron spectra from the PYTHIA calcu- 
lation with the high-pT-tuned parameters and modified 
Peterson function for charm quarks and beauty quarks 
are shown in the upper panel of Fig. [21 The PYTHIA 
calculation with the high-pr-tuned parameter set and 
modified heavy quark fragmentation function can simul- 
taneously describe the STAR direct open charm pT dis- 
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FIG. 4: Upper panels: Charged hadron sum pT distributions 
of D decay case and B decay case weighted by the fractions 
of electrons from D and B. The case that D meson decay 
contribution is always dominant is shown in the left panel. 
The case that B meson decay contribution become dominant 
at pt > 4.0 GeV/c range is shown in the right panel. Bottom 
panels: The ratios of pt > 5.0 GeV/c case topr > 2.0 GeV/c 
case are shown as solid circles. The ratios of pr > 3.5 GeV/c 
case to PT > 2.0 GeV/c case are shown as open circles. 



tribution and the STAR EMC electron data reported at 
the Quark Matter 2004 conference. Our calculation in- 
dicates that the contribution of electrons from B meson 
decays is not dominant for the measured px region up to 
8 GeV/c. 

The measurement of electron px distribution alone has 
a reduced sensitivity to the pr distribution of D mesons 
as shown in reference Within the large statistical 

and systematic errors of the electron data, it is possible to 
use high-pT-tuned PYTHIA calculation and the default 
Peterson function for both c and b quarks to yield an 
electron pT distribution similar to the measurement at 
high pt- In this shown in the bottom panel of 

Fig. 121 electrons from b quark decays are the dominant 
source for px greater than 4 GeV/c. 

In order to differentiate these two scenarios, we stud- 
ied the correlation between charm/beauty semileptonic 
decay electrons and charged hadrons in PYTHIA with 
these two sets of parameter in Peterson function. First, 
we select an electron decayed from B or D mesons at a 
given electron px cut, and make a cone around the direc- 
tion of the electron. Then we study the scalar summed 
Pt distributions of these charged hadrons in the cone 
(pt refers to the transverse momentum in the laboratory 
frame). Here the cone is defined by \rifi — Tyd < 0.35 and 
I'Ph ~ fel < 0.35 (77 is pseudorapidity and ip is azimuthal 
angle) . The summed pr distributions of charged hadrons 
at three different electron pr cuts are shown in Fig. |31 
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TABLE I: Fractions of B meson decay contribution at three 
electron pt cuts from two different PYTHIA calculations with 
high-pT-tuned parameters. Modified Peterson function and 
default Peterson function are each used for both charm quarks 
and beauty quarks. 



PT ( GeV/c ) 


2.0 


3.5 


5.0 


Modified Peterson function 


7% 


22% 


32% 


Default Peterson function 


17% 


54% 


69% 



The distributions are normalized to unity. The solid lines 
in Fig. O are electrons from B meson decays and the 
dashed lines are electrons from D meson decays. We 
found that the summed pT distributions do not change 
significantly between the default Peterson function and 
our modified fragmentation function. When electron pT 
cut goes from 2.0 GeV/c to 5.0 Gev/c, the summed pr 
distributions broaden for both B meson decay case and 
D meson decay case. However, the change is more drastic 
for the D meson decay case. 

The upper panel of Fig. 0] shows the sum of D decay 
case and B decay case weighted by fractions of electrons 
from D and B. Here the distributions are normalized 
such that the values in the first bin are the same. In the 
left upper panel, the fractions come from PYTHIA calcu- 
lations with the high-pT-tuned parameters and modified 
Peterson function for charm quarks and beauty quarks, 
while the fractions in the right upper panel of Fig. 0] are 
from PYTHIA calculations with the high-pr-tuned pa- 
rameters, but default Peterson function for charm quarks 
and beauty quarks. The values of the fractions are shown 
in Table I. When electron pT increases, if B meson de- 
cays become dominant, the summed pr distribution of 
charged hadrons changes slightly as shown in the right 
upper panel of Fig. 0] However, if the D meson decay 
contribution is always dominant, the summed pr distri- 
bution broadens drastically as shown in the left upper 
panel of Fig. ^ Ratios are shown in the bottom pan- 
els in Fig. ^ The solid cirles are the ratios of pt > 5.0 



GeV/c case over pT > 2.0 GeV/c case. The open cir- 
les are the ratios of pt > 3.5 GeV/c case over px > 2.0 
GeV/c case. Experimental observation of these ratios 
can be used to determine the relative significance of B 
and D decay electrons. 

In conclusion, we found that in order to match the 
STAR measurement oipx shape of D mesons from d-|-Au 
collisions at RHIC using PYTHIA Monte Carlo genera- 
tor, higher order QCD effects have to be taken into ac- 
count (PARP(64)=4) and a harder charm quark frag- 
mentation function must be used. The high-pr-tuned 
PYTHIA calculation with the modified fragmentation 
function can simultaneously describe the pr distribu- 
tions of D mesons and non-photonic electrons from semi- 
leptonic decays of heavy quarks. Because the measure- 
ment of electrons alone has a reduced sensitivity on the 
Pt shape of D mesons, the high-pr-tuned PYTHIA cal- 
culation with default Peterson fragmentation function 
for both c and b quarks can match the STAR electron 
measurement at high pT within large statistical and sys- 
tematic errors. This calculation predicts that the B de- 
cayed electrons dominate at the high pr {> 4: GeV/c) re- 
gion. However, in the proposed scheme of modified charm 
and beauty quark fragmentation, the beauty quark decay 
electrons are not a dominant contribution over the entire 
Pt region up to 8 GcV/c. Using PYTHIA as a guidance 
we found that the summed pt of charged hadrons within 
a narrow cone of heavy quark decayed electrons can be 
used to differentiate these scenarios. 
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